Quantum Billiards in Optical Lattices 
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We study finite two dimensional spin lattices with definite geometry (spin billiards) demonstrating 
the display of collective integrable or chaotic dynamics depending on their shape. We show that 
such systems can be quantum simulated by ultra-cold atoms in optical lattices and discuss how to 
identify their dynamical features in a realistic experimental setup. Possible applications are the 
simulation of quantum information tasks in mesoscopic devices. 



During the last decades, billiards have been the testbed 
of classical and quantum chaos [l|, due to their sim- 
plicity united to the richness of their displayed dynam- 
ics [3J. Theoretical evidences of the quantum mani- 
festation of chaos in billiards were first confirmed ex- 
perimentally in the spectral statistic of microwave res- 
onators 4j. quasi- two-dimensional superconducting res- 
onators 5(, and atom-optic billiards [6j. A very impor- 
tant step forward in the field occurred when it was re- 
alized that the properties of mesoscopic systems could 
be very sensitive, under appropriate circumstances, to 
the intestability properties of the underlying classical 
model [7|. We recall, as an example, the study of con- 
ductance fluctuations in quantum dots Q. More gen- 
erally, quantum billiards have shown to determine the 
dynamical properties of charge @, spin (lol | and entan- 
glement in nanostructures. The physics associated 
to quantum billiards has been shown very recently to be 
relevant in the study of graphene [IH ■ 

The ongoing interest in the study of quantum chaos 
stimulates the search for new physical systems where it 
is possible to experimentally study complex dynamical 
behaviour. In this Letter we propose to realize quantum 
billiards using optical lattices, which have been proved 
to be an excellent arena to study quantum many-body 
systems 
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Fundamental to our proposal is that optical lattices can 
operate as universal simulators 14j . i.e., by means of an 
appropriate dynamical control it is possible to reproduce 
the dynamics of any given spin Hamiltonian. Moreover, 
by means of the modelling of the form of the external trap 
it is possible to effectively define finite-size lattices. The 
class of billiards defined in this work are finite two dimen- 
sional optical lattices of given geometry, where collective 
excitations propagate and interfere as they back-reflect 
against the geometrical boundaries of the lattice. We 
therefore talk about spin billiards. This class shows a rich 
set of possible configurations, serving as a model system 
for different implementations: Depending on the system 
size, boundary conditions, lattice coordination number, 
and interaction Hamiltonian between the spins, one can 
either recover the known results on quantum billiards or 



model new physical systems showing original features. 
In the realm of cold atomic gases the distinction between 
regular and chaotic dynamics will appear in the momen- 
tum distribution of the atoms or in the fluorescence sig- 
nal. Present-day technology permits the simulation of 
the spin billiards introduced here. 

Optical lattices offer unique possibilities to simulate 
chaotic or integrable dynamics in a controlled way. The 
possibility to study billiards in this context gives a brand 
new perspective to a classic, and well studied, problem; 
numerous new questions can be addressed. On one side 
it is possible to explore the transition to chaos in a num- 
ber of different spin-Hamiltonian depending of its sym- 
metries. On the other side it is essential to understand 
the realization of the billiard, the measurement of rele- 
vant quantities and the sources of imperfections that may 
mask the physics we want to describe. We decide to first 
address this last topics. To this end we consider a model 
Hamiltonian which can be mapped on that of a particle 
hopping on a finite lattice (the billiard) . We will use the 
language of spin as it is natural in this case and it ap- 
plies also to those Hamiltonians where the mapping to a 
tight-binding model does not apply. 

The model - We consider a two dimensional 1/2-spin 
lattice with nearest-neighbor XX interaction in a trans- 
verse magnetic field. The Hamiltonian reads 



n = a 



<m,m / > 



v y 



(1) 



where a % a are the Pauli matrices, m = A4(i, j) is the com- 
posed index of the two dimensional qubits coding {i,j}, 
and the sum < to, to' > runs over nearest-neighbor spins 
on a square lattice with coordination number four (ex- 
cept at the boundaries) and free boundary conditions. 
We set A = 1 and h~ = 1. As TL commutes with the to- 
tal magnetization, we restrict to the subspacc with total 
magnetization equal to one (in this particular sector a 
mapping onto the single particle problem applies). The 
discrete space structure is a key- feature of optical lattices, 
we therefore need to re-discuss the effect of different bil- 
liard shapes on the dynamics. The first billiard under 
consideration is rectangular (billiard R, Fig.[IJ left). We 
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simulate the time evolution of a wave function initially 
peaked in an angle (single spin nipped, see Fig. [T] A). 
This situation can be realized by starting from a Mott 
insulator state with unit occupancy in a two-dimensional 
optical lattice [lij]. Two atomic hyperfine levels serve as 
the two pseudo-spin states, and all atoms can be pre- 
pared in one of the two by optical pumping. A laser can 
then be used to excite atoms from specific lattice sites to 
untrapped (continuum) states, removing them from the 
lattice. Its focus can be swept all along the border of 
a pre-defined region, leaving a regularly filled lattice of 
uniformly polarized atoms corresponding to the chosen 
billiard shape. At this point a Raman 7r-pulse on reso- 
nance with the transition between the two hyperfine lev- 
els can be used to flip the state of one atom at one of the 
billiard's vertices. In this particular geometry, shining 
the vertex atom with the edge of the laser spot will suf- 
fice, eliminating the need for subwavelength addressing. 
The simulation of the propagation of the resulting spin 
wave can then take place following a stroboscopic pro- 
cedure based on lattice-driven state-dependent collisions 
between neighboring atoms 14 , [lj| . Billiard S (Fig. [TJ 



right) is a quarter of a Bunimovich stadium and the ini- 
tial condition is, as before, localized at a boundary angle 
(see Fig. [T]B). In both cases, the initial condition reads 



|Vt=o) = |M(0,0)>. 



(2) 



Note that our analysis applies also to four- fold symmetric 
billiards, the symmetry axis of which passes through the 
site of the initial excitation A4(0,0). Our choice allows 
us to study the time evolution of an excitation neglecting 
the effects of central and axial symmetries. The initial 
condition ([2|) is such that the dynamics will be influ- 
enced also by very high energy levels, i.e., the dynamics 
is very far from being composed only by the low-lying 
excitations. On the contrary, in the continuum limit and 
starting with a different initial condition as, e.g., a Gaus- 
sian packet, one would recover the usual physics of elec- 
tronic billiards. Using the lattice on one hand leads to 
a discretization of the space, altering the geometric na- 
ture of curved edges (see Fig. QJ3, D, F); on the other 
hand, through the stroboscopic nature of the dynamic 
simulation it allows to "freeze" the system at a very well 
defined point in its time evolution for the purpose of 
state detection. The billiard dynamics is characterized 
by two time scales Tl and Ty. The first one is related 
to the characteristic length of the billiards L (~ 30 sites 
in our simulations), corresponding to the time needed 
for the first revival of excitations; the second timescale 
is given by the time needed to perform a swap between 
two neighboring spins, related to the inter-site coupling 
strength T\ = tt/ (4A). The relation of the two timescales 
is T L oc 2LT\. Fig. Q] C-F depict snapshots of the site 
excitation amplitude after time evolution at two differ- 
ent final times tf for billiards R and S: For tf > Tl/2 
the effect of different boundary shapes is already visible 




FIG. 1: Snapshots of the populations \ipt\ time evolution in 
the rectangular billiard (left) and stadium billiard (right), for 
t = 0+ (uppermost figures), t>T L /2 (middle) and t = 10 Tl 
(bottom). Color code goes from black (zero) to blue and red 
with increasing probability. 



(Fig. [T] C and D). For longer times, tf > T L the col- 
lective excitations spread all over the billiards, showing 
irregular profiles with no distinguishable features at first 
sight (Fig.[T]E and F). However, for tf cx n T L (n e N), a 
large revival at the initial site is still visible for the rect- 
angular billiard resulting from constructive interference, 
while this is no longer possible for the stadium billiard. 
These are the first signatures resembling chaotic and in- 
tegrable dynamics in billiards realized in optical lattices. 
In the following we demonstrate that this is indeed the 
case and that its characteristic features can be detected 
and quantified experimentally. 

We first check the level spacing statistics (LSS) for 
the two billiards R and S. Following Bohigas' conjec- 
ture we expect billiard R to show Poisson LSS, while 
billiard S should present something different due to the 
effect of level repulsion. Indeed, as shown in Fig. [21 we 
find that billiard R displays a well defined Poisson LSS 
(red squares). For billiard S, instead, we find a Semi- 
Poisson LSS typical of semi-integrable systems (also ap- 
pearing in the Anderson Metal- insulator transition) [l7| . 
The complete onset of chaos and the appearance of a 
Wigner-Dyson distribution is probably prevented due to 
the significant role still played by periodic orbits. A bet- 
ter convergence to the theoretical distribution can be 
obtained by considering defects, i.e., empty sites (see 
below). This allows a better statistics by averaging 
over Nr different configurations of defect probabilities 
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FIG. 2: Level spacing statistic for the billiard R (red squares) 
and billiard B (black circles). Dotted lines follow the Poisson 
(red) and Semi- Poisson statistics (black). Inset: LSS averaged 
over Nr = 10 different configurations of defects with Pd = 5 ■ 
1CT 3 (empty symbols) P D = 5-1CT 2 (full symbols). The black 
[red] dashed line follows the theoretical prediction P(s) = 
4sexp(-2s) [P(s) = exp(-s)]. 
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FIG. 3: Coarse grained Fidelity F3 (A,C) as a function of 
time for the integrable (black) and chaotic billiard (red) and 
the auto-correlation of the signals in the left figure (B,D) for 
P D = e = (upper panels) and P D = 5 ■ 10~ 3 ,e = 10~ 5 , 
Nr = 10 (lower panels). Insets: Magnification of the bigger 
figure. 



(Pd = 5 • 1CT 2 , 5 • 1CT 3 ), as shown in the inset of Fig. H 
Eventually, also billiard R displays Scmi-Poisson statis- 
tics due to the presence of defects (inset of Fig. [2J red 
full squares). 

The striking differences in the LSS discussed above are 
reflected in other features of the spin-billiard dynamics 
that can be measured experimentally as we show here- 
after. We consider in particular the momentum distri- 
bution and the fluorescence signal, as used for instance 
to detect single ions [18|. After its introduction by 



Peres [l9j, the survival probability or Fidelity F 

F=|(^=o|V t )| 2 , 



(3) 



has been very useful to characterize the transition to 
chaos (13 ■ Here the main issue, due to the lattice spacing 
coinciding with optical wavelengths, is single-atom spa- 
tial resolution. Therefore, we consider a Coarse Grained 
Fidelity (CGF) defined as 



<i+n,j+n> 

< E 

<hj> 



(4) 



i.e., the survival probability in a square region around 
the site {i, j}. The CGF can be obtained via fluores- 
cence measurements in optical lattices without single site 
addressing, which is at the edge of present day technol- 
ogy [HI- Even if the CGF fails in detecting the finest 
details of the dynamics it still captures the main differ- 
ences between the integrable and chaotic systems. In 
Fig. [3] A and C we show the CGF (n = 3) decay in the 
rectangular (black) and stadium (red) billiards for dif- 
ferent disorder settings together with their correspond- 
ing auto-correlation functions (Fig.[3]B and D). In spite 
of the random noise, both the CGF (Fig. [3] C) and its 



auto-correlation (Fig. [3] D) reveal the fundamental time 
scale Tl. Striking differences appear between R and S 
billiards: Periodic oscillations persist up to times of the 
order of tf Tl in the integrable case (black), while 
in the chaotic case (red) a damping shows up on a time 
scale t >Tl revealing a rapid decay of correlations. 

Finally, we investigate the momentum distribution of 
the magnetization in the two billiards at a final time t f : 
The results are reported in Fig. [4] A and B for the rect- 
angular (left) and stadium billiards (right). Again the 
results show striking differences: In the integrable case 
the number of frequencies relevant to the wave function 
are much less than in the chaotic case. Moreover in the 
former case there is a structure (even though quite com- 
plex) in the spectrum that is absent in the latter one. 

Experimental Implementation - As mentioned in the 
introduction spin billiards can be studied experimentally 
in optical lattices following the idea of the Universal 
Quantum Simulator 14]. The measurements of the CGF 
and momentum distribution can be performed via fluo- 
rescence and timc-of- flight methods respectively [l|l [H| . 
The signatures of chaos we arc interested in arise from 
the evolution of a single-spin excitation corresponding to 
single-particle signals. In order to attain a sufficient atom 
number resolution to detect it, an average over different 
realizations is required. Under ideal conditions, since the 
evolution is fully deterministic, this should not be a prob- 
lem and the result would be fully reproducible. However, 
in a real experiment, errors might introduce differences 
between repetitions which may result in a corrupted out- 
put. The major error sources are two: (i) imperfections 
in the realization of a square lattice with hundreds of sites 
with uniform occupation number one, and (ii) side effects 
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FIG. 4: 2D FT of the Wave function F(u) x ,u y ) for the rect- 
angular (left) and stadium (right) billiards at time tf = 10 Tl 
for P D = e = (A,B); P D = 5 • 10~ 3 ,e = W~ S ,N R = 10 
(C,D); and P D = 10 _2 ,e = lO" 5 ,^ = 10 (E,F). Color code 
is the same of Fig. [T] 



of the parabolic magnetic field trapping the atoms in the 
studied region of the optical lattice [22j . We can model 
these errors via the presence of defects or "holes" in the 
spin billiard (missing atom in a optical lattice site) and 
errors in the gates performed to simulate the dynamics. 
The Hamiltonian ^ is then replaced by 

Hi = \ J2 (<< n +<<')+£ ^x*+j>r, (5) 

<m,m'> m 

where m = M.H(i,j) takes into accounts the presence 
of defects and e(t) fluctuates in [0, e] with flat distribu- 
tion. We repeat the previous analysis accounting for ex- 
perimental errors with typical values Pr> = 5 • 10~ 3 and 
e = 10 -5 averaging over Nr different configurations [22| . 
Fig. [3] (lower panels) shows the coarse grained fidelity F 3 
and the auto-correlation as a function of time R and S 
billiards in presence of experimental errors. In this case, 
it is more difficult to distinguish between integrable and 
chaotic dynamics; however, a careful inspection can still 
reveal differences. As before, in the integrable case the 
auto-correlation revivals last for longer times tf 3> Tl 
and their visibility is greater than in the chaotic case. 
More clear signatures are found again in the momentum 
distribution (Fig. [4j lower panels): The structures in the 
frequency domain lasts for very long times in the inte- 
grable billiard (even if slightly blurred) while they disap- 
pear in the chaotic case. 

Finally we would like to highlight the possible devel- 
opments along the lines presented here: The study and 



simulation of weak localization, quantum hall effect, dis- 
order effects, quantum information protocols, entangle- 
ment dynamics, and the role of different Hamiltonian 
and/or parameter regimes. We also point out that a 
similar analysis could be performed for alternative exper- 
imental setups as, e.g., lattices of coupled cavities [23j . 
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